Quadratic Form Expansions for Unitaries 



Niel de Beaudrap^, Vincent Danos^, Elham Kashefi'^, Martin Roetteler^ 



2 



^ IQC, University of Waterloo 



School of Informatics, University of Edinburgh 
^ Laboratoire d'Informatique de Grenoble 
NEC Laboratories America, Inc. 



Abstract. We introduce techniques to analyze unitary operations in terms of 
quadratic form expansions, a form similar to a sum over paths in the computa- 
tional basis when the phase contributed by each path is described by a quadratic 
form over R. We show how to relate such a form to an entangled resource akin to 
that of the one-way measurement model of quantum computing. Using this, we 
describe various conditions under which it is possible to efficiently implement a 
unitary operation U, either when provided a quadratic form expansion for U as 
input, or by finding a quadratic form expansion for U from other input data. 



1 Introduction 

In the one-way measurement model [1*21, quantum states are transformed using 
single-qubit measurements on an entangled state, which is prepared from an 
input state by performing controlled- Z operations on pairs of qubits, including 
the input system and ancillas prepared in the |+) state. This model lends itself 
to ways of analyzing quantum computation which do not naturally arise in the 
circuit model, e.g. with respect to depth complexity [3 | and discrete structures 
underlying unitary operations 16181 . In this article, we present another result of 
this variety, by introducing quadratic form expansions. 

Definition 1. Let V be a set ofn elements, and 1,0 '^V be (possibly intersect- 
ing) subsets. For a binary string x G {0, 1}^, let x/ and xo be the restriction 
ofx to those bit-positions indexed by elements of I and O, respectively. Then a 
quadratic form expansion is a matrix-valued expression of the form 

U = ^ Yl ^''^^"^ |xo)(x/|, (1) 
xe{o,i}^ 

U : 0.2^ — > Ti®^, where Q is a real-valued quadratic form on x, and C £ C 

Quadratic form expansions bear a formal similarity to a representation of a prop- 
agator of a quantum system in terms of a sum over paths. For a unitary U given 
as in ([T]), the outer product |xo)(x/| essentially specifies a particular coeffi- 
cient, in the row indexed by the substring x/ and the column indexed by xq: 
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the amplitude of the transition between these standard basis states is propor- 
tional to a sum of complex units specified by x/, xo, and the auxiliary variables 

V eV \{I u O). 

Representations of unitary transformations as sums over paths is a well- 
developed subject in theoretical physics (see e.g. II4I5I ): and a representation 
of unitaries as a sum over paths was used in [9l to provide a simple proof of 
BQP C PPQ However, there are also examples of quadratic form expansions 
which arise without explicitly seeking to represent unitaries in terms of path 
integrals: the quantum Fourier Transform over 1^2" can readily be expressed 
in such a form, and quadratic form expansions for Clifford group operations 
are implicit in the work of Dehaene and de Moor lITTl . as we will describe in 
Section O 

Given such an expression for a unitary U, we show how to obtain a decom- 
position of U in terms of operations similar to those used in the one-way mea- 
surement model. Using this connection, we demonstrate techniques involving 
quadratic form expansions to efficiently implement a unitary operator, when the 
coefficients of the quadratic form satisfies certain constraints related to "general- 
ized flows" (or gflows) [8] or Clifford group operations. In particular, we exhibit 
an 0(n^/ log n) algorithm to obtain a reduced measurement pattern (an algo- 
rithm in the one-way model) for Clifford group operations from a description of 
how they transform the Pauli group, based on the results of ifTTl . 

2 Connection to the one-way model 
2.1 Review of the one-way model 

We can formulate the one-way measurement model as a way of transforming 
quantum states in the following way. Given a state ona set of qubits / (the 
input system), we embed / in a larger system V, where the qubits of F \ / are 
prepared in the |+) oc |0) + |1) state. We then perform entangling operations on 
the qubits of V, by performing controlled- Z (denoted f^Z) operations on some 
sets of pairs of qubits. (These operations are symmetric and commute with each 
other, and so we may characterize the entangling stage by a simple graph G 
whose vertices are the qubits of V: we call this the entanglement graph of the 
procedure.) We then measure each of the qubits of V in some sequence, except 
for some set of qubits O C y (the output subsystem) which will support a final 
quantum state. We may represent the measurement result for each qubit v by a 

' Unitaries were expressed in |9| in terms of pathis whose phase contributions are described by 
cubic polynomials over Z2; comments made in Section VI of that paper essentially anticipate 
quadratic form expansions with discretized coefficients. We describe how their techniques 
provide a means of constructing quadratic form expansions from circuits in Appendix lAl 
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bit E {0, 1} which indexes the orthonormal basis states of the measurement. 
The measurement basis for each qubit may depend on the results of previous 
measurements, but without loss of generality may be expressed in terms of a 
"default" basis which is used when all preceding measurements yield the result 
0. Depending on the measurement results, a final Pauli operator may be applied 
to the qubits in the output subsystem oH 

In the original formulation of the one-way measurement model, the mea- 
surement bases were described by some axis of the Bloch sphere lying on the 
XY plane, which is sufficient for universal quantum computation. It is also easy 
to prove that restricting this to states which are an angle 6* G f Z from the X axis 
is sufficient for approximately universal quantum computation [12|. While it is 
reasonable to extend beyond this for choices of measurement bases [7 ], we wiU 
only need to consider measurement bases from the XY plane. 



2.2 Phase map decompositions from quadratic form expansions 

Consider a unitary U given by a quadratic form expansion as in ([Hi, where the 
quadratic form Q is given by 

Q(X) : 



(2) 



for some angles {Ouv}u vev > ^^'^ where the sum includes terms for u = v. Note 
that (5(x) can be expressed as an expectation value (x| H |x), where H is a 
2-local diagonal operator: 



H 



E 

{u,v}CV 



Then we may decompose U as follows: 



U oc |xo)(x| e'-^|x)(x/| 



xe{o,i}^ 



E 



|yo)(y| 



.ye{o,i}^ 



iH 



X/ 



.xe{o,i}^ 



oc Roe'" Pi 



(4) 



where Pi is a unitary embedding which introduces fresh ancillas (indexed by 
V ^ P = V \ I) initialized to the | +) state, and Rq is a map projecting onto 
the |-|-) state for all qubits mO^ = V \ (tracing those qubits out afterwards). 



" The reason for using the same variables V, I, and O for these sets of (labels for) qubits as for 
the sets in Definition[T]will become apparent in the next section. 
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Equation (01) is a phase map decomposition flO] for U: that is, it expresses 
U in terms of a process of postselecting observables, as follows. Decompose H 
into terms Hq, Hi, and H2, where Hq consists of the 1-local terms on the qubits 
of O, Hi consists of the 1-local term on the remaining qubits, and H2 contains 
the remaining terms from We then have U oc Rq e'^'^e^^^e^^^ Pj. Note 
that e^^ci and e*^^ are simply single-qubit Z rotations applied to the elements 
of O and respectively, where in each case the qubits v in those sets are ro- 
tated by an angle 9^^ . Then, the composite map Rq = Roe^^^ projects each 
the state of each qubit v ^ onto the vector |0) + e"*^"" |1) for each v G O*^. 
We then have U = e^^o p^^ 

is a decomposition of U into the 
prepai^ation of some number of |+) states, followed by a diagonal unitary op- 
erator consisting of two-qubit (fractional) controlled-Z operations, followed by 
post-selection of states on the Bloch equator for v G O^, and (unconditionally 
applied) single-qubit Z rotations on the remaining qubits. If 9uv G {0, vr} for all 
distinct u,v £ V and for u = v £ O, the above describes precisely the action of 
a measurement-based computation in which the qubits u G O"^ are measured in 
the eigenbases of observables of the form M{—6yy) = cos(0^,^,)X — sm{6yy)Y, 
in the special case where all measurements result in the +1 eigenstate (which 
we may label with the bit S„ = 0). 

If we are able to extend the above into a complete measurement algorithm, 
with defined behavior when not all measurements yield a specific outcome, we 
obtain a measurement-based algorithm for U : we discuss this problem in the 
next section. Conversely, from every measurement based algorithm, we may 
obtain a quadratic form expansion: 

Theorem 1. Every unitary operator on n qubits may be expressed by a quadratic 
form expansion with \I\ = \0\ = n, and where the quadratic form has coeffi- 
cients 6uv G {0,7:} for all cross-terms XuX^ and — vr < 9^.^ ^ i: for all terms 
x1 . Furthermore, any unitary can be approximated to arbitrary precision by 
such an expansion where we further require 9^ S ^Z. 

Proof. From | H] (and using the notation of that article), the measurement pat- 
tern M~'^ EuvNy performs the unitary transformation J {a) = [ } ] 
for a G M, from the state space of a qubit u to that of a "fresh" qubit v. These 
operations generate SU(2), and generate a group dense in SU(2) if we restrict 
toaG fZ.bv Ifni. 

For any n qubit unitary U, there exists a measurement pattern composed of 
such patterns together with two-qubit controlled- Z operations (which we denote 
AZ) which implements U. Let G be the entanglement graph of this pattern, and 
I and O be the qubits defining the input space and output space (respectively) of 
the measurement pattern. By [6], in this measurement pattern, the probability of 
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every measurement resulting in the +1 eigenvalue (i.e. S,; = for all v G O^) 
is non-zero. Then, U oc Rq Pj , where 



H = - 

uv£E{G) 



|l)(l|n®|l>(l| 



By dll), this yields a quadratic form expansion for U, with 

(5(x) = ^ TTXuXv - ^ 



(5) 



(6) 



For a quadratic form expansion approximating U, it is sufficient to consider 
measurement patterns approximating U using angles a„ G f Z. □ 



2.3 Measurement Pattern Interpolation 

As we remarked above, the connection from quadratic form expansions to phase 
map decompositions may allow us to obtain an implementation for U, provided 
we can determine how to adapt measurements in case the measurements for 
qubits f G O"^ do not all yield the result S„ = 0. 

In a measurement pattern performing N measurements, the computation 
may follow any of 2^ branches, corresponding to the different combinations of 
measurement results. Let us call the branch in which every measurement pro- 
duces the result = the positive branch of the measurement patternJl With- 
out loss of generality, we may restrict our attention to patterns where no clas- 
sical feed-forward is required in the positive branch: then, the positive branch 
of a measurement pattern is characterized by the geometry {G, I, O) of the pat- 
tern (where G is the entanglement graph of the measurement algorithm, and 
I,0 y{G) are the sets of qubits defining the input/output space of the pat- 
tern), and the angles a = {a„}^g(^c defining the measurements to be performed. 

To extend the description of the positive branch of a measurement algorithm 
into a complete measurement algorithm performing a unitary is the subject of 
the following problem: 

Measurement Pattern Interpolation (MPI). For input data {G, I, O, a), de- 
scribing a unitary embedding U as the positive branch of a measurement pattern 
with geometry (G, I, O) and performing measurements a, determine if there a 
measurement pattern *P with geometry {G, /, O) which performs the transfor- 
mation U. 

^ This choice of terminology refers to all measurements yielding the +1 eigenvalues of their 
respective observables Af(— 
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This problem is open, and seems to be difficult in general. We may attempt to 
make the problem easier by considering a more restricted problem: 

Generic Measurement Pattern Interpolation (GMPI). For an input geome- 
try {G,I,0), determine if there exist measurement patterns ^P(a) parameter- 
ized by a choice a of measurement angles, each with geometry (G, /, O), such 
that the pattern ^(a) performs a unitary embedding for all a. 

GMPI addresses, in an angle-independent manner, the subject of the structure of 
measurement patterns which perform unitary transformations. A special case of 
the GMPI which has been solved are those geometries (G, /, O) which have a 
"generalized flow" (or gflow), which are the ''yes" instances of GMPI such that 
the patterns ^(a) yield maximally random outcomes on all of their measure- 
ments [8|. The following is the definition of gflows in I.13il . for measurements 
restricted to the XY plane Q 

Definition 2. Given a geometry (G, /, O) for a measurement pattern, a gflow 
is a pair (g, =<;), where g is a function from O'^ to subsets ofl'^ and ^ is a partial 
order, such that the following conditions hold for all u and v in the graph G: 



where odd(S') is the set of vertices adjacent to an odd number of elements of S. 

Here, u ^ v essentially represents, for two qubits u and v, that v is measured 
no earlier than u; a gflow then specifies an ordering in which the qubits are to 
be measured (with the function g providing a description of how to adapt later 
measurements). Mhalla and Perdrix |T3 1 present an algorithm which determines 
if a geometry has a gflow in this sense in polynomial time, which in turn yields 
a polynomial time solution to the GMPI for that case. As a result, any instance 
of the MPI where the geometry (G, /, O) has a gflow can be efficiently solved. 

A different special case of the Measurement Pattern Interpolation problem 
which has been solved is that where the measurement angles ai^e restricted to 
multiples of | (or slightly more generally, where the measurement observables 
are Pauli operations). In this case, as noted in Q, no measurement adaptations 
are necessary, and the corrections can be determined via the stabilizer formal- 
ism na. 

In the following section, we apply these solutions to special cases of the 
MPI to describe how to synthesize implementations for a unitary U given by a 
quadratic form expansion. 

* The original definition of gflows in {8) also allows for YZ plane and XZ plane measurements, 
which do not play a role either in our analysis or in \l3i . 




(7a) 
(7b) 
(7c) 




u £ odd{g{u)) 
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3 Synthesis via measurement pattern interpolation 

In order to apply the partial solutions to the MPI described above, it will be 
useful to define the following: 

Definition 3. For a quadratic form expansion 

1 ^ ^iQiy^) |xo)(x/| where Q(x) = ^ OuvXuXv , (8) 
xe{o,i}^ {■^M'^v 

the geometry induced by the quadratic form is a triple {G, I, O) , where G is 
a weighted graph with vertex-set V, edge-set {uv \ u ^ v and Ouv 7^ 0} , and 
edge-weights Wg{uv) = ^™/vr. 

Because we can require — vr < Ouv ^ for all n, v G V , we may without 
loss of generality restrict G to have edge- weights — 1 < Wg{uv) ^ 1. We will 
assume that this holds for the remainder of the article, and speak of edges being 
either of unit weight or fractional weight. 

In this section, we consider the problem of synthesizing an efficient imple- 
mentation of unitaries U in terms of the geometry induced by a quadratic form 
expansion for U by reduction to the solved cases of the Measurement Pattern 
Interpolation problem discussed in the previous section. 



3.1 Measurement pattern synthesis via gfiows 

Consider a geometry (G, /, O) induced by a quadratic form expansion for a 
unitary embedding U, where G has only edges of unit weight: then (G, /, O) is 
also a geometry for a measurement pattern. To obtain a measurement pattern for 
U, it suffices to find a gflow for (G, /, O): in that case, by Theorem 2 of [8], for 
any choice of measurement angles a = {a^j^g^c, we may consider the pattern 



n 



1 ii£odd(g(«)) 




veg{u) 



\{Nu 



(9) 



where the left-hand product may be ordered right-to-left in any linear extension 
of the order ^, and ~ denotes the adjacency relation of G. This pattern thus 
steers the reduced state after every measurement to the state which would occur 
if the result had been the +1 eigenvalue. Every branch of the pattern then per- 
forms the same operation as the positive branch, and so the pattern implements 
a unitary operation U. To obtain a pattern in standard form (with corrections 
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only on output qubits), it is sufficient to propagate the corrections to the left, 
absorbing them into the measurement bases. 

In lfT3l . an 0{n^) algorithm is provided to determine whether or not a ge- 
ometry (G, /, O) has a gflow where every qubit is to be measured in the XY 
plane (and obtain one in the case that one exists), where n = \ V{G)\. The mea- 
surement pattern of ^ can be constructed in time O(n^) by first producing a 
pattern where corrections undo byproduct operations after each measurement, 
commuting these corrections to the end, and simplifying; the resulting pattern 
will have 0(n) operations each with complexity 0(n). Thus: 

Theorem 2. For a unitary embedding U given as a quadratic form expansion 
with geometry {G, /, O) with unit edge-weights, there is an 0{n^) algorithm 
which either determines that (G, I, O) has no gflow, or constructs a measure- 
ment pattern consisting ofO{v?) operation^ implementing U (using measure- 
ment angles of arbitrary precision), where n = \ V{G)\. 

3.2 Circuit synthesis via flows 

A geometry (G, /, O) which has fractional edges lies, at first glance, outside of 
the domain of the Measurement Interpolation Problems described above. How- 
ever, given a quadratic form expansion with such a geometry, we may still be 
able to synthesize a circuit for a unitary U represented by that expansion by 
considering flows, which con^espond to gflows where the function g maps each 
vertex u E to a singleton set: we may say (/, ^) is a flow if and only if 
(5/, ^) is a gflow, where gf{v) = {f{v)]. 

Geometries which have flows are a solvable special case of the GMPI, where 
the resulting measurement patterns are very "circuit-like". Specifically, the pos- 
itive branch of a measurement pattern whose geometry has a flow can be repre- 
sented by a circuit with the following characteristics [6|: 

- edges of the form v f{v) for v ^ correspond to J(— a^,) gates on some 
wire, separating two wire segments which we label v and f{v); 

- edges uv G E{G) for u ^ f{v) and v / f{u) correspond to f\Z operations 
acting on the wire segments labelled by u and v; 

- wires whose initial segments are labelled by vertices of / accept arbitrary 
input states, while those labelled by vertices P \ img(/) take input |+). 

In the above formulation, the edges of the form v f{v) can be interpreted as 
implementing single-qubit teleportation, in which case a fully entangling unitary 

^ These operations may involve measurement angles of arbitrary precision. A corresponding ap- 
proximate measurement pattern may use 0{n^ +npolylog(n/e)) operations by the Solovay- 
Kitaev Theorem WAi . where e is the precision of the coefficients of Q. 
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is important in order to transfer the information of the "source" qubit to the 
"target" qubit upon measurement. However, considering the analysis of [6|, it 
is not important that the edges of the second kind above be fully entangling 
operations: using such edges to represent fractional powers of i\Z will also yield 
unitary circuits. This motivates the following definition: 

Definition 4. Suppose {G, I, O) is a geometry of a quadratic form expansion 
for a unitary transformation U. We may say that (/, is a fractional-edge 
flow /or (G, /, O) if it is a flow for that geometry, and for all ah G E{G) with 
Wciab) < I, we have f{a) / b and f{h) / a. 

If (G, /, O) has a fractional-edge flow, we may synthesize a circuit from a 
quadratic form expansion for U using the description above, where edges ab 
of fractional weight correspond to AZ'^'^^"*) gates on the wire segments la- 
belled by a and b rather than simple gates. We will make use the following 
easily verified Lemma to consider how to compose/decompose quadratic form 
expansions: 

Lemma 1. Let Ui, U2 be matrices given by quadratic form expansions 



(/) IfVi r\V2 = h = Oi, then U2U1 = ^ E e*QiW + *Q2(x) |xo2)(x/- 



(//) IfVi and V2 are disjoint, then U10U2 = e^'^i W + *Q2(x) |xo)(x7| , 

where I = h LI I2 and O = Oi U O2 . 

We prove the circuit construction given by inducting on the number of edges 
of fractional weight. For the base case, if (G, /, O) has no fractional-weight 
edges at all, we may synthesize a circuit for U as above, as it corresponds to a 
normal measurement pattern with a flow, and so falls under the analysis of JH. 
We may then induct for geometries with fractional edge-weights if we can show 
we can decompose the geometry into ones with fewer fractional edge-weights. 

For any arbitrary fractional edge ab G E{G) and each each z G O, we may 
define m{ab, z) to be the maximal vertex v G V{G) in the ordering =^ subject to 
z being in the orbit of v under / (that is, z = f^{v) for some i ^ 0), such that 
at least one ofv ^. aorv ^b holds. For a set 5 C V{G), let G[S] represent the 
subgraph of G induced by S (i.e. by deleting all vertices in G not in S). Then, 
define the following subgraphs of G, and corresponding geometries: 




(10) 



In the following, G = C1C2, and sums are over {0,1} 



VxVJV2 



X 
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(G^^^O) {Gl,I,V2) {G2,V2,V2) {G3,V2,0) 




V3 O O 

Fig. 1. Illustration of the decomposion of a quadratic form expansion about an 
edge ab, expressed in terms of geometries. V2 is a set of maximal vertices under 
the constraint of being bounded from above, by the vertices a and b, in a partial 
order ^ associated with a fractional-edge flow. Arrows represent the action of 
the corresponding fractional-edge flow function, /. 

- Let V2 be the set of vertices m{ab, z) for each z e O^: it is easy to show 
that a, b € V2. Let G2 = G[V2], and let ^2 = (G2, F2, V2). 

- Let Vi be the set of vertices u G V{G) such that u ^ v for some v ^¥2; let 
Gi = G[Vi] \{uv\u,ve V2] ; and let Gi = {Gi,I, V2). 

- Let V3 be the set of vertices u E V{G) such that u ^ t> for some v V2; let 
G3 = G[V3]\{uv\u,veV2}; and let = (G3, V2,0). 

This decomposes the geometry (G, /, O) into three geometries with fractional- 
edge flows, as illustrated in Figure [U 

Let Qi be a quadratic form on {0, 1}^^ consisting of the terms XuX^ of Q for 
u £ Vi or V ^ Vi, but not both; Q2 be a quadratic form on {0, 1}^^ consisting 
of the terms XuX-a of Q for distinct u,v £ V2; and similarly let be defined on 
{0, 1}^^, and consist of the remaining terms of Q. Then Qi, Q2, and define 
quadratic form expansions for some operations Ui, U2, and U3 (respectively) 
with geometries Gi,Q2, and (respectively). 

- U2 in particular will be a product of operations AZ^^^'^'"^ for distinct u,v G 
V2 , as it is a quadratic form expansion whose input and output indices co- 
incide. Then U2 can be represented as a circuit with a wire for each u E V2, 
with fractional controUed-Z gates AZ^'^C"'') for each edge uv G E{G). 

- Both Qi and Q-^ have fractional-edge flows, but fewer fractional edges than 
(G, /, O). By induction, Ui and C/3 are also unitary embeddings, and have 
circuits with wire-segments connected by J {9^) gates (where 9^ are the co- 
efficients of the terms x1 in each quadratic form) and possibly fractional f\Z 
gates (as in the case for U2)- 

- In the circuits described above, the terminal wire-segments for Ui and (a 
subset of) the initial wire-segments for C/3 have the same labels as the wires 
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Fig. 2. The geometry for the quadratic form expansion of the QFT for Z32, and 
the con^esponding circuit due to (iT]. In the geometry (on the left), input vertices 
are labelled by circles, output vertices by lozenges, and fractional edges are 
labelled with their edge-weights. 

for C/2 ■ The composite circuit for U?jU2Ui can then use these labels to arrive 
at a unified labelling of its' wire-segments. 

Because <5i(x^^ ) + Q2{^V2 ) + Q^i'^v-^ ) = Q(^) for all x G {0, 1}^ by con- 
struction, the composite operation U3U2U1 can differ from U by at most a scalar 
factor by Lemma [TJ so the circuit obtained implements the operation U. 

In |[T3l . an 0{kn) algorithm is provided to determine whether or not a ge- 
ometry (G, /, O) has a flow, and obtain one if it exists, where n = \V{G)\ and 
k = \0\. For each edge uv, we may check whether one of Wg{uv) = 1 or 
[n / f{v) and v 7^ f{u)\ holds: if all edges satisfy this constraint, the circuit 
described above is well-defined. By iterating through the vertices of V{G) in an 
arbitrary linear extension of ^ , we may construct the circuit described above 
can be constructed in time 0{m), and the size of the resulting circuit will also 
be 0{m), where m = \E{G)\. By an extremal result |15|, any geometry with a 
flow has m ^ kn: thus, the total running time of this algorithm is 0{kn). 

In the case |/| = |0|, a flow function / is unique if it exists, by ll20l : so in 
this case, if (G, /, O) has a flow but there is an edge v f{v) of fractional weight, 
there is no fractional- weight flow for (G, /, O). We then have: 

Theorem 3. For a unitary transformation U given as a quadratic form expan- 
sion with geometry (G, /, O), there is an 0{kn) algorithm which either deter- 
mines that (G, /, O) has no fractional-edge flow, or constructs a circuit consist- 
ing of 0{kn) operation^ implementing U, where n = \ V{G)\ and k = \0\. 



^ These operations may consist of J(a) gates and fractional /\Z gates of arbitrary precision. A 
corresponding circuit using a finite elementary gate set may be of size 0{kn polylog(fcn/£)) 
by the Solovay-Kitaev Theorem 614i . where e is the precision of the coefficients of Q. 
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Example. The Fourier Transform over is given by the matrix formula 



^ n — 1 71 — 1 



'2 

x,ye{0,l} 



(11) 



which is a quadratic form expansion; its quadratic form can be given by 

n— 1 n—l—h 

2"-! 



Q(x,y) = 2^ 2^ -i^;^T^Xhyj. (12) 

h=0 j=0 



This has a fractional-edge flow for all n. Figure |2] illustrates this geometry for 
n = 5, and the circuit (due to ETl ) which may be synthesized from it. 

3.3 Synthesizing measurement patterns for the Clifford group 

If a quadratic form expansion has a geometry whose edges all have unit weight, 
and its' other coefficients are multiples of |, then it corresponds to the positive 
branch of a measurement pattern which measures only X or y observables. A 
measurement pattern of this sort, if it performs a unitary operation, performs a 
Clifford group operation in particular. 

An algorithm of Aaronson and Gottesman |fT9l can produce a circuit of size 
0{v? / \ogn) in classical deterministic time 0(n^/ log n) for a Clifford group 
operation U acting on n qubits, from a description of how U transforms Pauli 
operators by conjugation. By converting the circuit into a measurement-based 
algorithm, and performing the graph transformations of [ 1 8 1 to remove auxiliary 
qubits, we may obtain a pattern of at most 3n qubit^ in time 0(n^/ log n). 
Building on the results of ifTTl . we show how to classically compute such a 
minimal pattern in time 0{n^/ log n) by solving the MPl for a quadratic form 
expansion for U. 



Obtaining a quadratic form expansion. For the sake of completeness, we 
outline the relevant results of ifTTl . Define the following notation for bit-flip and 
phase-flip operators on a qubit t out of a collection {1, . . . , n}: 

Pt = , Pn+t = Zt . (13) 

Let diag(Af) G represent the vector of the diagonal elements of any square 
boolean matrix M; and let d(M) = diag {M^ [ ° \"] M) G Z|" for a 2n x 2n 

' In I?), Clifford operations on n qubits are described as having minimal patterns for are de- 
scribed as requiring at most 2n qubits; however, this only holds up to local Clifford operations 
on the output qubits. 
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matrix M over Z2. Then, we may represent an n qubit unitary [/ by a 2n x 2n 
boolean matrix C and a vector h e {0, 1}^", whose coefficients are jointly 
given by 



UPtU^ 



(14) 



for each 1 ^ t ^ 2n. (Note that the factor of is only necessary to ensure 
that the image of Pt is Hermitian, and does not serve as a constraint on the value 
of C as a matrix.) We will call an ordered pair (C, h) a Leuven tableau for a 
Clifford group element U if it satisfies (fT4l) j^l 

Provided a Leuven tableau (C, h) for a Clifford group operation U, iflTl pro- 
vides a matrix formula for U which we may obtain for U, as follows. Decom- 
pose C as a block matrix C = [ ^ |J ] with n x n blocks, and then find 
invertible matrices Ri , R2 over Z2 such that k^^GR2 = [ 1 J for 

some 

r < n (using e.g. the decomposition algorithm of |22| to obtain Ri and R2 in 
terms of elementary row operations). Then, define the matrices 



Ell E12 
E21 E22 



Ri ER2 , Ri 



i?i , R2 



E. 
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1,. 



R. 



2 5 



(15) 



where En is taken to be a block of size (n — r) x (n — r). We may then obtain 
the block matrices 



In-r £^12 -Pll F12 

E2I E22 F21 F22 

Hii H12 

1^ H21 H22 



Rj 
Rj"' 



C 



r; 



and use these to construct the n x n boolean matrices 



(16) 



Mbr = 

Next, define 



Fii + E12H21 E12 
EJ2 E22 



Mbc 



H^i 
H21 H22 



(17) 



= diag(Mfer) , dbc = diag(Mbc) , 

Lbr = lower (Mfcr + dft^d^) , L^c = lower {Mbc + dbcdj^) 



(18) 



' Note that the block matrix [ hi is similar to a destabilizer tableau as defined in 1191 . 
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where lower (M) is the strictly lower-triangular part of a square matrix M (with 
all other coefficients set to 0). Finally, define il^ = [[] i^] and U:^ = 1„ — 77^ 
for the sake of brevity, and lej^ 



[ 1„ ] h + diag ( [7^2 ^77^] Lb,, [7?2 ^77^ 



[0 7?2 



h + fl,"' diag 7?, 



Lbc + HrMbc 



+ (n^+ HrMbc Lbr + MbcH, 



R2 



(19) 
(20) 



Then Theorem 6 of [17] states that the unitary operation U for the Clifford 
operation characterized by (C, h) is given by the matrix formula 



U 



1 



X6e{0,i}"-'- 



be be J ^ 



xc,x,e{0,l}'- (_i)«x,. + d^c^i.c) \R-^Xbr){R2^Xbc + t\ 



(21) 



where ^br = [ x^ ] and ^bc = [ xc ] are n bit boolean vectors. 

The formula in (|2T]) shows strong similarities to a quadratic form expansion. 
In particular, consider disjoint sets of indices Vb, Vr, and Vc, with \ Vb\ = n — r 

and \Vr\ = \Vc\ = r. Let V = Vb U Vc U Vr, I = Vb U Vc, and O = Vb U Vr, 



and define the following notation for x G {0, 1}^ : 

G {0, 1}^ , xo 



X/ 



Xf, 













Xb 










xw_ 



G{0,1}° , 

Q(x) = TT(xJ)LbrXo + X^IIrXi + xj Lbc:>Cl + x/hfech^c^/ 



vr 



- ( Xodbrd5,.xo + xj dfecd^^x/ 



(22) 

(23) 



Then, (|2TI ) is equivalent to 
1 



/2r 



X6{0,1}^ 



(24) 



which is essentially a quadratic form expansion sandwiched between two net- 
works of controlled-not and X gates. To obtain a simple quadratic form expan- 
sion, we would like to perform a change of variables on x/ and xo; but this 
cannot be done as 7 and O intersect at Vb, and the changes of variables do not 



' The vector formulas given here for t and hi,c rnay be obtained by repeated application of 
Theorem 2 of 1171 . 



Vt = i n o 
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• — 


• 









/ O' 

Fig. 3. Illustration of geometries arising from quadratic form expansions yield- 
ing the same matrix. On the left is a geometry whose inputs and output inter- 
sect; on the right is a geometry from an equivalent quadratic form expansion, 
constructed so that the input and output indices are disjoint. 

necessarily respect the partitioning of I and O with respect to this intersection. 
However, we may add auxiliary variables in order to produce an expansion with 
disjoint input and output indices. Note that 

I2 = 5x„.Jx2)(xi| = \ (-l)"^"^+"^"Ma^2)(xi| (25) 

xe{o,i}^ xe{o,i}^ 

where 5x^y is the Kronecker delta. Let Va and Vy be disjoint copies of Vh, and 
set V' = V y^Va^^ Vy and O' = Vy U Vr- Writing and x;,' for the restriction 
ofxG {0, l}^'toV;andT4/, we then define 

G {0,1}^ , xo' 



X6 



Xfc' 



G{0,1} 



O' 



(26) 



Q'(X7, Xa, Xo') = 7r(x5,L6,.Xo' + xJ/TIrX/ + xjLfecX/ + h^^X/ 



+ TTX; 



U- n-r 





Xa + vrx, 



O' 







x„ 



- ^(df[,xo' + d1,i^i^ . (27) 

Note that the difference between the expressions for Q' and Q is essentially that 
all instances of xo have been replaced with xq/ (which is independent from x/), 
and the presence of the terms involving x^. (This manipulation is illustrated in 
Figure [3] as a transformation of geometries.) We therefore have 

e^«W|i?ixo>(i?2'x, + t| 

xe{o,i}^ 

= <5x„x,,e^^'(''^'0'''o')|i?ixo'>(i?2ix, + t| 

= ^ E e*'^'('^^'^-^o')|^^x^,)<^-ix, + t| . (28) 
xeio,!}"" 
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Substituting the final expression of (1281 ) into (|24l) and performing the appropriate 
change of variables, we have 

[/ = ^ ^ g^Q'(K.(x,+t),x„,i?,-Xo,)|xo,)(x,| . (29) 

xe{o,i}^' 

Note that the quadratic form of the expansion in ( [29l) has only angles O^v 
which are multiples of |, with 9u^ G {0, tt} for u v. This then represents the 
positive branch of a one-way measurement pattern on the geometry (G', /, O') 
of the quadratic form expansion of |29j using only X or y basis measurements, 
and having only n — r auxiliary vertices. 



Interpolating the measurement pattern. We can augment this to a measure- 
ment pattern by applying the techniques of the stabilizer formalism |[T6l to the 
stabilizer code generated by the operators K{v) = W^^w v £ P 

(where again ~ is the adjacency relation of G), as follows. To obtain the final 
correction, we do classical pre-processing simulating the evolution of the state 
space when we perform one measurement at a time. For each measured qubit u, 
there is an associated correction which we may perform immediately after 
the measurement if we obtain the result = 1. We store for each qubit v two 
boolean formulas (3^ and jy, representing the X and Z components of the ac- 
cumulated corrections to be performed on v. When v is measured, the pending 
X corrections will affect the result of any Y measurement, and the pending Z 
corrections will affect the result of any X or Y measurement, in each case by 
exchanging the significance of the two measurement outcomes^ Just prior to 
the (simulated) measurement of 6^, let 5„ = jyifv is to be measured with an X 
observable, and 6^ = 0^+ Iv if v is to be measured with a Y observable. Thus, 
upon measuring v, the following operations are accumulated into the corrections 
which must be performed: 

- For every qubit w where acts with an X or y operation, we must add 

+ 6v to (3u,; 

- For every qubit w where acts with a Y or Z operation, we must add 

Sy + 6y to Jw. 

This accounts for the accumulated corrections due to the measurement of v and 
every preceding measurement which affects it. By simulating measurement for 
all of the qubits in in this way, we obtain boolean formulae for the correc- 
tions on O in terms of the results of the measurements: the correction to be 



This can be described in terms of signal shifting, as described in yjj. 
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performed for some w S O is Z'^" , for Pw and 'jw constructed after all of 
the (simulated) measurements. To obtain {3^ and 7^^, for all w £ O in this way 
takes time O(n^). 

It is easy to show that the resulting measurement pattern is irreducible by the 
techniques of fTSl, by the following argument. Let A denote the set of auxiliary 
vertices corresponding to the bit positions of x^: note that in the measurement 
pattern, these are all to be measured with the observable X, and are adjacent 
only to the input/output variables x/ and xo0 To eliminate a vertex u G A 
using the methods of [ 1 8 ] on the geometry induced by the quadratic form ex- 
pansion, we would have to identify an output variable 60 G O adjacent to x, 
and apply the graph transformation in |[T8l Proposition 1]. This would result in 
a geometry where has the former neighbors of u in G (and in particular is 
not adjacent to any more removable vertices), and where a local Clifford (which 
is not a Pauli operator) must be applied to bo after the entangling procedure. 
Because b^ is not adjacent to any other auxiliary qubit after this transforma- 
tion, the local Clifford cannot be undone or made into a Pauli operator by e.g. 
another vertex removal; then, except by extending the computational model to 
allow for corrections which are local Clifford operations, performing the local 
Clifford can only be done by introducing an auxiliary qubit (or rather, a new 
output qubit following bo, making the latter an auxiliary qubit). Thus: 

Theorem 4. For an n-qubit Clifford group operation U given in the form of a 
Leuven tableau, there is an 0(n^/ log n) algorithm which produces a minimal 
one-way measurement pattern for U. 

The ability to obtain a quadratic form expansion representing a reduced 
measurement pattern yields a more efficient algorithm to find totally reduced 
Clifford patterns than from using existing techniques to obtain one via the cir- 
cuit model. The quadratic form of ( [29l ) can be found from a Leuven tableau 
(C, h) in time 0{n^ / logn), which is dominated by the time required to com- 
pute i?i and R2- To contrast, an approximately optimal quantum circuit for a 
Clifford group operation {i.e. consisting of 0(n^/logn) gates) can be found 
from a Leuven tableau in time 0{rv' / log n) by transforming it into a destabi- 
Uzer tableau, and then applying the algorithm of 1 19 ]. To obtain a measurement 
pattern from such a circuit by composing the patterns for each gate, removing 
vertices opportunistically (with each removal taking time 0{v?)), requires time 
0{n^ I log n). Thus, making use of quadratic form expansions provides us with 
a faster algorithm to obtain reduced measurement patterns for Clifford group 
operations. 

" There are no square terms x\ for n G A or cross-term x^Xt, for u,v ^ A before the change 
of variables in J28t), and the change of variables itself does not introduce any. 
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4 Conclusions and Open Problems 

We have introduced quadratic form expansions, and developed techniques which 
suggest that they may be useful for synthesizing efficient implementations for 
unitary operations. We described conditions under which implementations may 
be efficiently found for unitaries specified by quadratic form expansions; and 
we showed how quadratic form expansions leads to more efficient algorithms 
for obtaining reduced patterns for Clifford operations in the one way measure- 
ment model. 

In the introduction, we mentioned that quadratic form expansions are sim- 
ilar in form to a sum-over-paths representation of unitary operations, which 
is a well-developed subject in theoretical physics. This raises the question of 
whether the techniques developed here are useful e.g. for developing algorithms 
to simulate physical systems. It is not known whether the solved cases of the 
Measurement Pattern Interpolation problem correspond to natural (in the more 
Uteral sense) unitaries expressed as sums over paths: this question, and how to 
extend the solved cases of the MPI to include propagators for interesting physi- 
cal systems, remain open. 
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A Quadratic form expansions as sums over paths 

Let [G, I, O) be the geometry of a quadratic form expansion, as defined on 
page |7] In tiie special case when {G, I, O) has a fractional-edge flow as de- 
fined in Section IT2I the quadratic form expansion corresponds exactly to a sum 
over paths as described in [9|, for the elementary gate set of H, Z*, and a^*, 
where t ^ R {i.e. admitting arbitrary Z rotations and fractional controlled- Z 
gates). In order to demonstrate the sense in which quadratic form expansions 
are sums over paths in this case, and because it represents a reasonably simple 
algorithm for converting quantum circuits into quadratic form expansions, we 
now present an alternate proof of Theorem [T] based on the techniques of ||9l. 
That any quadratic form expansion with geometry with a fractional-edge flow 
can be constructed in this way follows by reversing the construction below. 

Proof of Theorem [TJ Consider a quantum circuit implementing U exactly, 
using the operations H, aZ*, and Z*. Enumerate the wires of the circuit from 1 
to k, and for each wire 1 ^ j ^ k, introduce a path label xj for the input end 
of the wire, corresponding to an input bit Xj G {0, 1}. We set / = {1, . . . , A:}. 
Divide each wire into segments, bounded on each end by either a Hadamard 
gate, the input terminal of the wire, or the output terminal. We label the wire 
segments with path variables: for the segments at the inputs, we apply the labels 
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Xj for j S /, and we introduce new path variables to label the remaining wire 
segments. Computational paths in the circuit are then described by setting all 
of the the path variables xi - ■ ■ Xn collectively to some particular binary string 
in {0, 1}". The phase contribution of each paths, governing how they interfere 
to produce an output state for any given input state, is described by a function 
(/9(x) depending the gates of the circuit as follows: 

(/) For every Hadamard gate on a single wire, with a path variable x^ labelling 
the segment preceding the Hadamard and a path variable Xj labelling the 
segment following the Hadamard, we add a term x^xj . 

(//) For every a^* operation between two wires, with a path variable Xh labelling 
the segment of one wire and Xj labelling the segment of the other wire in 
which the aZ* operation is performed, we add a term tx^Xj . 

{Hi) For every operation on a wire segment labelled with a path variable Xj , 
we add a term (Because the path variable Xj ranges over {0, 1}, the 
extra power of 2 has no effect.) 

In particular, the function (/^(x) is a quadratic form, where without loss of gen- 
erality the coefficients may be constrained to — 1 < t ^ 1. The phase of a given 
path, described by a bit-string x G {0, 1}", is then given by (-l)'^(^) = e''^'^^''). 
Each path also has an associated amplitude of 2^'''/'^, where r = n — A; is the 
number of Hadamard gates in the circuitF^ 

Let O be the set of indices j such that some wire is labelled by the path- 
variable Xj at its' output end. Then, the initial points of computational paths 
are described by bit-vectors a G {0, 1}^ , and the terminal points of paths are 
described by b G {0, l}*^ . The coefficients C/b,a can then be given as the sum 
of the contributions of all paths beginning at x/ = a and ending at xo = b : 

xe{o,i}" 

x/=a 
xo=b 

which is an expression of the coefficients of C/ as a quadratic form expansion. 

To obtain a proof of Theorem [TJ it is sufficient to note that without loss of 
generality we may restrict ourselves to using A^* gates only for t = 1 to imple- 
ment U exactly; and that to implement U to arbitrary precision, it suffices to use 
Z* gates where t is restricted to multiples of |. □ 



Although it is quite reasonable to consider to be simply a polynomial over R, in terms of the 
descriptions used in Section VI of |9|, one may consider <^ to be a polynomial over the ring 
R/2Z. If we restrict to t G jZ, we may simplify this to the finite ring Zg by multiplying all 
of the coefficients by 4, and using it to describe powers of \/i rather than of —1. 



